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Q I Abstract 

For equidistant discretizations of fractional Brownian motion (fBm), the prob- 
■ abilities of ordinal patterns of order d = 2 are monotonically related to the Hurst 

p ^ i parameter H. By plugging the sample relative frequency of those patterns indicat- 

ing changes between up and down into the monotonic relation to H, one obtains the 
Zero Crossing (ZC) estimator of the Hurst parameter which has found considerable 
attention in mathematical and applied research. 
i | In this paper, we generally discuss the estimation of ordinal pattern probabilities 

in fBm. As it turns out, according to the sufficiency principle, for ordinal patterns 
^ , of order d = 2 any reasonable estimator is an affine functional of the sample relative 

00 I frequency of changes. We establish strong consistency of the estimators and show 

' them to be asymptotically normal iov H < |. Further, we derive confidence intervals 

. for the Hurst parameter. Simulation studies show that the ZC estimator has larger 

' I variance but less bias than the HEAF estimator of the Hurst parameter. 

o ; 

00 
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^ ' 1 Introduction 

Probabilities of ordinal patterns in equidistant discretizations of fractional Brownian mo- 
tion (fBm) have first been analyzed by Bandt and Shiha Ordinal patterns represent 
the rank order of successive equidistant values of a discrete time series. For example, for 
three successive values there are six different possible outcomes of the rank order, which 
we call ordinal patterns of order d = 2. 

As Bandt and Shiha have shown, for equidistant discretizations of fBm the distribution 
of ordinal patterns is stationary and does not depend on the particular sampling interval 
length. Further, the probabilities of ordinal patterns of order d = 2 are all monotoni- 
cally related to the Hurst parameter, with the probability of ordinal patterns indicating 
changes from up to down and from down to up, respectively, being strictly monotonically 
decreasing in H. 

The estimator of the Hurst parameter obtained by plugging the sample relative frequency 
of changes between up and down into the monotonic functional relation to H has been 
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known for some time, running under the label 'Zero Crossing' (ZC) estimator since changes 
between up and down correspond to zero crossings of the first order differences. 

Note that more generally, but not focussing on fBm, Kedem [11] has considered the 
estimation of parts of the autocorrelation structure of a stationary Gaussian process by 
counting (higher order) zero crossings, that is, zero crossings of the first (or higher) order 
differences. 

The ZC estimator is asymptotically normally distributed for < |, a result essentially 
due to Ho and Sun [10], who proved that the sample relative frequency of changes is 
asymptotically normally distributed in this case. CoeurjoUy [6] has resumed properties of 
the ZC estimator including strong consistency. In the sequel, the applicability of the ZC 
estimator has been examined by Markovic and Koch [17] and Shi et al. [20], including 
simulation studies as well as the application to hydrological and meteorological data. 

The distribution of ordinal patterns is the base of ordinal time series analysis, being a 
new fast, robust and flexible approach to the investigation of large and complex time 
series (see Bandt [2] and Keller et al. [T3l [12] ) . From the viewpoint of ordinal time series 
analysis, the estimation of ordinal pattern probabilities in stochastic processes is of special 
interest . 

This paper is structured as follows: Sec. [2] is devoted to a general discussion of the 
estimation of ordinal pattern probabilities in fBm where we consider estimators based on 
counting the occurence of ordinal patterns in realizations. Obviously, the sample relative 
frequency of ordinal patterns is an unbiased estimate of the probability of the patterns. As 
Theorem m shows, averaging the sample relative frequencies of ordinal patterns and of their 
'time' and 'spatial' reversals, one obtains estimates which are strictly more concentrated 
in convex order. Hence, we restrict our following considerations to such 'reasonable' 
estimators. Notice that Theorem H] has important consequences for the estimation of 
functionals of ordinal pattern distributions such as the permutation entropy (see [2], |12j). 

For ordinal patterns of order d = 2 one obtains two different estimators which can both 
be expressed as affine functionals of the sample relative frequency of changes. By Lemma 
El we establish strong consistency and asymptotical unbiasedness for the estimators of 
bounded continuous functionals of ordinal pattern probabilities. Theorem [8] states that 
estimators of ordinal pattern probabilities are asymptotically normal for if < |. 

In Sec. [31 we consider the estimation of the probability of a change by the sample rel- 
ative frequency of changes in some detail. We give a formula for the precise numerical 
evaluation of the variance of the sample relative frequency of changes as well as asymp- 
totically equivalent expressions. Based on these results, confidence intervals for the Hurst 
parameter are provided in Sec. [H 

In Sec. [5l we compare the previous results to findings for simulations of fBm. As it turns 
out, the confidence intervals obtained from the ZC estimator in Sec. [H cover the true 
unkown value of the Hurst parameter for about 95 per cent of the cases, even for small 
sample sizes and for the values of the Hurst parameter larger than |, where asymptotic 
normality of the ZC estimator is not necessarily expected to hold. Compared to the HEAF 
estimator which estimates the Hurst parameter by plugging the sample autocovariance 
into a monotonic functional relation to H, the ZC estimator has larger variance but much 
less bias, in particular for small sample sizes. 

Notice that the results given in this work for the increments of fBm similarly apply to 
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FARIMA(0,(i,0) processes. Furthermore, the statistical properties of the estimates of ordi- 
nal pattern probabilities do not only apply to fBm, but also to montonic transformations 
of fBm. In particular, the estimates of ordinal pattern probabilities are invariant with 
respect to (unknown) non-linear monotonic transformations of processes. 



2 Estimating ordinal pattern probabilities 

2.1 General aspects 

Let {Q, A) be a measurable space equipped with a family (P^)^ge of probability measures 
with G non-empty. The subscript indicates that a corresponding quantity (E^, Var^, 
etc.) is taken with respect to F^. For the special case of fractional Brownian motion 
where = ]0, 1], we write H instead of -d. 

Let (Xfc)jtgNg be a given real-valued stochastic process defined on {Q,A). Define the 
process of increments (Ifc)fcgN by Y^. := X^. — Xk_i for /c G N. For (i G N let Sd denote the 
set of permutations of {0, 1, . . . , li}. 

Definition 1. For c? G N let the mapping tt : R'^^^ Sd be defined by 

7r((xo,Xi,...,Xrf)) = ( ° I I ■■■ ^ ) =: (ro, n, . . . , r^) 

y To ri r2 ... Td J 

for (xo, Xi, . . . , Xd) G W'-^'^ with the permutation (ro, ri, . . . , r^) of {0, 1, ... , d\ satisfying 
Xd-ro > Xd-ri > ■■■> Xd-ra, and r;_i > ri if Xd-n.^ = Xd-r, for I = 1,2, ... ,d. For /c G Nq 
we call 

Ud{k) := 7i{{Xk,Xk+i,...,Xk+d)) 
the (random) ordinal pattern of order d at time k. □ 

The permutation vr((xo, Xi, . . . , Xd)) describes the rank order of the values Xq,Xi, . . . , Xd, 
where in case Xd-i = Xd-j for < i < j < d the 'earlier' Xd-j is ranked higher than Xd-i- 

Notice that if (Xjt)jtgi^g is pairwise distinct, that is, Xi ^ Xj P^-a.s. for all i,j G No 
with i j and G 0, then, apart from sets with probability zero for all G 0, Ild{k) 
generates the same cr-algebra as the rank vector {Rq, Ri, . . . , Rd) of (X^, Xk+i, . . . , Xk+d) 
given by Rj = Yfi=Q l{Xfc+,>Xfc+a for j = 0, 1, . . . , (see [I5], p. 286). 



Stationarity. For -i? G let = denote equality in distribution with respect to P^. A 
stochastic process {Zk)k£T defined on {Q,A) for T = N or T = No := N U {0} is called 
stationary iff with respect to each ■(9 G 

{Zki, Zk2, ■ ■ ■ , Zk„) — [Z^-^^i, Zk2+i, ■ ■ ■ , Zk^+i) 
for all ki, k2, . . . , kn & T with n G N and for all / G N. 

In fact, Ild{k) only depends on (Y^+i, ^A:+2, • • • , ^fc+d) for all (i G N and k G Nq. In 
particular, let 

7f((yi,2/2, • • • ,Z/d)) = (ro,ri, . . . ,rd) 



3 



be the unique permutation of {0, 1, ... , d} for 1/2, ■ ■ ■ , yd) G M'^ such that 

d-ro d—ri d—r^-i d-r^ 

j=i j=i j=i j=i 

and ri_i > ri if Yl'j=i'^ yj = Ei=? yj for / = 1, 2, . . . , d. Obviously, 

7r{{xo,Xi,...,Xd)) = TT{{xi-Xo,X2-Xi,...,Xd-Xd-l)) 

for all (xo, Xi,...,Xd) e M'^+\ and hence Ild{k) = vf ((F^+i, • • • , Yk+d)) for all /c G Nq. 
This immediately yields the following statement. 

Corollary 2. If(Yk)keN is stationary then (lldik)) ^^^^^ is stationary. 

Space and time symmetry. For d G N let the mappings a, (3 from 5*^ onto Sd be 
defined by 

a{r) := {rd,rd-i,...,ro), p{r) := {d - Tq, d - r^, . . . , d - Vd) (2) 

for r = (ro, ri, . . . , rd) G Sd- Geometrically, a(r) and /5(r) can be seen as the spatial and 
time reversal of r, respectively (see Figured]). Let the set r be defined by 

r := {r,a{r), P{r), P o a{r)} (3) 

with o denoting the usual composition of functions. As a o /?(r) = jjo a{r) and a oa{r) = 
(3 o f3{r) = r, one has «(r) = f3{f) = f. Clearly, if s G r for r, s G Sd, then s = f. This 
provides a division of each Sd into classes, consisting of 2 or 4 elements. For d = 1 the 
only class is Sd = {(0, 1), (1,0)}, for d = 2 there are the two classes {(0, 1,2), (2, 1,0)} 
and {(0, 2, 1), (2, 0, 1), (1, 2, 0), (1, 0, 2)}, and for c/ = 3 there are 8 classes. Note that for 
d > 3 classes of both 2 and 4 elements are possible. 




r = (0,2,1) a(r) = (1,2,0) /3(r) = (2, 0, 1) /3 o a(r) = (1, 0, 2) 

Figure 1: The ordinal pattern r = (0, 2, 1), its spatial reversal a(r), its time reversal /3(r), 
and its spatial and time reversal (3 o a{r). 

For d E N and n G N let (Sd)"' denote the n-fold Cartesian product of Sd- Define the 
mappings A, B from (Sd)^ onto (5*^)" by 

A((r«,r(2),...,r("))) := (a(r«), ^(r^^)), . . . , a(r("))) , 
B((r«,r(2),...,r("))) := (/5(r(")), /3(r("-^)), . . . , /3(r«)) 

for {r^^\ r^'^\ - - - , r^^^) G (Sd)""- Further, let id denote the identity map on (5'^)". Note 
that {id. A, B, B o A} together with o forms an Abelian group. 
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We say that (yfe)feeN is symmetric in space and time iff with respect to each -i? G 

and = (n„,n„_„...,nj 

for all ki, ^2, • • • , fcn e N with n G N. 

Lemma 3. // {Xk)k&no ^■^ pairwise distinct and {Yk)km is symmetric in space and time, 
then the distribution o/II = (11^(0), 11^(1), . . . , ^di'^ — 1)) is invariant under the group of 
transformations {{id, A, B, B o A}, o| on (5"^)" with respect to for each -i? G 0, i.e. 

n ^ A(n) ^ B(n) ^ BoA(n). 

Proof. If (Xfc)fegN,) is pairwise distinct then, according to the definition of tt, 

Y2,...,Yd)),..., niiVn, . . . , Y^+n-i))) 
= A{{ni{-Y^, -Y,, -Y,)), . . . , ^-Y^, -Y^+^, -Y,+n-i)))) 
= B{{n{iYd+n-i, . . . , Y^+u . . . , ni{Yd, . . . , F2, Y,)))) 
= B o A((?F((-Frf+„_i, . . . , . . . , ^{{-Y,, -Y,, -Y^)))) 

P^-a.s. with respect to each -i? G (see ([1])). Further, symmetry in space and time 

of (Yk)km implies that (Fi, F2, • • • , {-Yi, -Y2, -Ya+n-i), • • • , >2, ^1) 

and (— Yrf+„_i, . . . , — 12, — ^i) have the same distribution with respect to for each -i? G 0, 

and hence the statement follows. □ 



A Rao-BIackwellization. For fixed r ^ Sd with d G N consider the functional Pr{-) 
defined by Pr{^) '■= ^^(11^(0) = r) for G 0. Obviously, if (yfe)A;6N is stationary then the 
statistic 

n— 1 

P^,n = P^J'^) ■= -$^l{n.(fc)=r} (4) 

fc=0 

of n = (11^(0), 11^(1), . . . , Ild{n — 1)) is an unbiased estimate of Pr{-), that is E^(p^^) = 
Pr{'^) for all G and n G N. If, additionally, {Xk)kmo i^ pairwise distinct and (Yk)k£n 
is symmetric in space and time then, according to Lemma El the statistic 

p^,^ = p^ju) := ^(^p^jn)+p;:jA{u))+p;jB{u))+p^jBoAin))') (5) 

of n is a Rao-Blackwellization of (see Theorem 3.2.1 in [18j)- This proves the following 
Theorem. 

Theorem 4. Let r G IJdeN'^'i- V {^k)ke^o pairwise distinct and {Yk)k&i is stationary 
and symmetric in space and time then the estimate of pr{-) is unbiased and more 
concentrated in convex order than p^,^, that is 

E^((^(p^^,p,(t9))) < E^(^(p;r„,p,(^))) (6) 

for all ^ E Q with respect to each function ip : [0, 1] x [0, 1] [0, oo[ with (p{p,p) = and 
^p{-,p) being convex for every p G [0, 1]. 
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According to the strictness of the Jensen inequahty for strictly convex functions, if 
(Xfc)fcgNo is pairwise distinct and {Yk)k£¥i is stationary and symmetric in space and time, 
one has strict inequahty in ([6]) whenever (p{-,p) is strictly convex for every p G [0, 1] and, 
additionally, P^(p^^ 7^ pjT^) > (see p/S], Theorem 3.2.1). Since (■— p)^ is strictly convex 
for every p G [0, 1], in particular 

Var^(p;;^) < Vartf(p;r^) 

in this case. Note that for r E Sd with d = 1, the estimator simply estimates the 
constant functional Pr{-) = |. 



Ergodicity. Consider the measurable space {Q', A') := (R^, i3(M^)) of infinite sequences 
of real numbers and let the mapping T : Q' ^ Q' he defined by T{uj') = {uj2,uj'^, . . .) for 
u' = {uj[,LJ2, . . .) G Q'. The process (yA:)^^^ is called ergodic iff (yA:)^^^ is stationary and, 
additionally, for every A E A' such that T~^{A) = A one has P^((lfc)fcGN G A) = or 
Ptf((>^fc)fceN e A) = 1 for each G 6. 

As the next Lemma shows, if {Yk)km is ergodic then the estimators of continuous and 
bounded functionals of ordinal pattern probabilities are strongly consistent and asymp- 
totically unbiased. 

Lemma 5. Let r G IJrfeN ^d- If (Xk)keN is ergodic and h : [0, 1] ^ M is continuous then 

Jirn/i(p^^) = h[pr{'&)) 

P^-a.s. for all d eQ. If h is continuous and hounded, then with respect to each d E Q 

\immJh{p;^S) = KpM))- 



Proof. For fixed -(9 G 9 let the probability measure /x on {Q! , A!) = (M^, i3(M^)) be defined 
by n{A) := P^((Yfc)fcgN e A) for A E A'. Further, let the mapping T : Q' ^ fl' he given 
as above and define / : fi' ^ M by 



1 for IT {{uj[,uj'2,. . .,uj'a)) E r 
else 



for u' = {uj[,uj2, . . .) G Q'. Obviously, / is Borel-measurable and J^, \ f\dfi < 00, hence, 
according to Birkhoff's Theorem, if {Yk)km is ergodic then 



^ n— 1 

hmp^^ = lim - V/(r^((n)fceN)) = E^(/((n)fceM)) = Pr{^) 



j=0 



P^-a.s. (see [7], Theorem 1.2.1) and the first statement follows since h is continuous. If 
additionally h is bounded then 

Jim E^(/i(p^j) = E^(^ lim h{p;j^ = h{pr{d)) 
according to the dominated convergence Theorem. □ 
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2.2 Specialization to fBm 



In the following we specialize the considerations of Subsection 12.11 to equidistant dis- 
cretizations of fractional Brownian motion (fBm). We start with the definition of fBm. 

Definition 6. For H G ]0, 1] let (i?()f:))fg[o,oo[ be a mean-zero Gaussian process on a 
probability space (17, A, P) with the covariance function 

Cov(5(t),5(s)) = ^(t'^ + s'^- Var(S(l)) 

for s,t G [0,oo[. Then (i?(t))tg[o,oo[ is called fractional Brownian motion (fBm) with the 
Hurst parameter H. □ 

It is well-known that each fBm possesses a modification with P-a.s. continuous paths 
(see [9]). We are only interested in the distribution of fBm, so let for the rest of the 
paper {^h)h&\q,i] be a family of probability measures on a measurable space (f2, ^), and 
(-B(t))(g[o,oo[ be a family of real-valued random variables on {VL^A) such that (i?(t))(g[o,oo[ 
measured with respect to ¥h is fBm with the Hurst parameter H. E.g., (i?(t))ig[o,oo[ can 
be defined as the identity on the set of continuous functions on [0, oo[. 

For the sampling interval length 5 > consider the equidistant discretization 

{Xl)k&io '■= {B{k6))keNo 

of fBm. Further, let := - Xl_^ for k e N. Note that fBm is if-self-similar (see 
[9]), that is, for every H e]0,1] and for all a G ]0, oo[ it holds 

(5(at))te[o,oo[ ^— {(i^B(t))teio,ool 

with ^= denoting equality of all finite-dimensional distributions with respect to Fh here. 
Since it holds tt ((yi, 2/2, • • • , Ud)) = ((a'^Z/i, a^2/2, • • • , a^Vd)) for aW a > 0, H e ]0, 1] and 
{yi,y2, ■■■,yd) e R'^, one has 

(7r((>^fe+l, • • • , yk+d)))kmo (7f((>^fc+l' • • • ' ^fcVd)))fe6No 

with respect to each H G ]0, 1] for all 6 > 0. This means that the distribution of ordinal 
patterns for equidistant discretizations of fBm does not depend on the particular sampling 
interval length. Therefore, in the following we only consider the discretization {Xk)keNo = 
{XDkeNo with the increment process {Yk)km = (Xl - Xl_-^)keN- 

Note that, similarly to the sampling interval length, one can show that the particular scal- 
ing Var//(i?(l)) of fBm has no effect on the distribution of ordinal patterns for equidistant 
discretizations. Hence, we always assume the case of standard fBm where VaiHiBll)) = 1 
for all H G ]0,1]. 

According to Definition [6l VainiXi — Xj) > for all i,j G No with i ^ j and H G ]0, 1], 
and since Xj — Xj is Gaussian with respect to Fh for all H G ]0, 1], the stochastic process 
{Xk)keNo is pairwise distinct. Further, obviously (Yfe)A:eN is mean-zero Gaussian with 
respect to Fh for each H G ]0, 1], and by Definition [6] one obtains 

PH{k) := CoYH{Y^,Y^+k) = + - 2|fc|2^ + |fc - (7) 
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for all H e ]0, 1] and G Nq, that is, the stochastic process (yfc)fcgN is stationary. Fur- 
thermore, for all n G N and ki, k2, ■ ■ ■ , kn G N the random vectors (Yfe^, Y^j, . . . , Y^^), 
(— Yfci, — Yfe2, . . . , — Yfc„) and (Yfc„, . . . , Ifca; ^fcJ have the same covariance structure, hence 
{Yk)k^f>i is space and time symmetric. Consequently, the conclusion of Theorem H] applies 
to the estimation of ordinal pattern probabilities for equidistant discretizations of fBm. 

As the following Lemma shows, one has strict inequality in (P) whenever ip{-,p) is strictly 
convex for every p G [0, 1] and if G ]0, 1[ (compare to the remark after Theorem Hj). 

Lemma 7. If H e ]0, 1[ then P//(p^„ 7^ JV^„) > for all r G IJ^^pj Sd and n eN. 

Proof. Let r E Sd with G N and n G N be fixed. It is easy to construct a sequence of per- 
mutations (^r^^\r^'^\ . . . ,r(")) G (Sd)"^ such that the (n + rf — l)-dimensional Lebesgue mea- 
sure of HLi { (2/i>?/2, • • • ,?/n+d-i) e I 7r((?/fe,2/fc+i, • • .,yk+d-i)) = r^'')} is strictly 
positive and 

tlr ■ tt{A;G {0, l,...,n-l}|r('=) = r} ^ ^ {'^ ^ {0, 1, . . . , n - 1} | r^'^) G r}. 

According to ([7]), if if G ]0, 1[ then {Yi, Y2, . . . , Yn+d-i) is non-degenerate Gaussian with 
respect to and, consequently, Fh{{Yi,Y2, . . . ,Yn+d-i) G A) > for each Borel-set 
A C R'^^"^^ with strictly positive {n + d — l)-dimensional Lebesgue measure. Since 
~ ff'l^serPs^n (^^^ ® ^^'^ ®) Statement follows. □ 

Note that if if = 1 then P/f ^ p^,J > if and only if r G {(0, 1, . . . , rf), (rf, . . . , 1, 0)}. 

It is well-known that for each H G ]0, 1[ the spectral distribution function of {Yk)km is 
absolutely continuous (see [1]) which is a sufficient condition for (Yfc)^^^ to be ergodic 
with respect to Pjy (see [7], Theorem 14.2.1). In the case H = 1 one has = kXi 
P/^-a.s. for every k eNq (see [9]), hence 



ifre[j,^^{Sd\{{OA,...,d),{d,...,l,0)}) 
\ ifrGaeN{(0,l,---,^^),(f^,---,l,0)} 



Pjy-a.s. for all n G N. Consequently, the conclusions of Lemma [5] apply to the estimation 
of ordinal pattern probabilities for equidistant discretizations of fBm. 



Asymptotic normality. Next, we discuss asymptotic normality of the estimators of 
ordinal pattern probabilities for equidistant discretizations of fBm. We restrict our con- 

si deration to the 'reasonable' estimators as provided by Theorem |H Let — ^ denote con- 
vergence in distribution with respect to P/f, and A^(/i, a^) be the normal distribution with 
mean G M and variance cx^ > 0. Given a mean-zero Gaussian vector Z = (Zi, Z2, . . . , Zd) 
on a probability space (f2',^',P), and a Borel-measurable function / : M*^ — »• M such that 
Var(/(Z)) < cx), define the rank of / with respect to Z by 

rank(/) := min |k G N : There exists a real polynomial g : M*^ ^ M 

of degree k with E([/(Z) - E(/(Z))] g(Z)) ^ O}, 

where the minimum of the empty set is infinity. We write f{k) ~ g{k) for mappings f,g 
from No onto M and say that / is asymptotically equivalent to g iff limfc^oo f{k)/g{k) = 1 
where ^ := 1. Note that for puik) as defined in (JTj), one has the asymptotic equivalence 

PH{k) ~ i7(2i7 - l)A;2^-2 (9) 
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for each H G ]0, 1] (see [H]). Further, we write f{k) = 0{g{k)) iff sup^,gj^j^ f{k)/g{k) < oo, 
and fik) = o{g{k)) iff Wuik^^ f {k) / g{k) = 0. 

Theorem 8. If H < \ then 

{Ye.THip?J)-'-{p;^^-p,.{H)) ^ N{0,1) 

for all r e Sd with d 

Proof Let r e Sd with G N \ {1} be fixed and define / : M'^ ^ M by 
f[[yr,y2,...,yd)) - | q ^^^^ 

for (yi, 1/2, • • • , e M'^. Write Y = (Fi, F2, • • • , We first show that rank(/) > 1 
with respect to Y and ¥h for all H G ]0,1]. Obviously, / is Borel-measurable, and 
Var//(/(Y)) < 00 for all H G ]0, 1]. Now, let z G {1, 2, . . . , be fixed. Because {Yk)km is 
symmetric in space and time, the random vectors (Fi, 1^2, • • • , Yd) and (— Fi, —Y2, . . . , —Yd) 
have the same distribution with respect to P// for all H G ]0, 1], and hence, according to 
the definition of a (see (E])) 

E.//(l{?F(Y)=s}'i^i) + E/f (l{^F(Y)=Q{s)}Fj) 

= Eh(1{^(y)=s}V;) -E^f(l|^(_Y)=s}(-l"i)) = (10) 

for each s G 5*^. Since Yi is mean-zero Gaussian with respect to ¥h for all H G ]0, 1], one 
has Ej/(/(Y)) W^uiYi) = 0, and this together with ( ITOl) yields in case fjf = 2 

EH([/(Y)-EH(/(Y))]r,) = EH{l{^iY)=r}Yi) +EH{l{^Y)=a(r)}Y^ = 

with respect to each H G ]0, 1], and in case {jr = 4 

E^,([/(Y)-EH(/(Y))]r,) = EH{lMY)=r}Yi)+EH{l{:,(^Y)=air)}Yi) 

+ ^H{'^{w{Y)=l3(r)}Yi) +EH{l{n{Y)=l3oa(r)}Yi) =0 

with respect to each H G ]0, 1]. Consequently, rank(/) > 1 with respect to Y and for 
all H G ]0, 1]. According to (ED, \pH{k)\"'''^^^'^ = 0{k^^~^) for each H G ]0, 1], thus 

oo 

^|Ph(A:)|™'^(^^ < oo 

fc=0 

for if < |. Since, by definition one has \ Yl=l fiiYk+i, • • • , Yk+d)) = Pt,^ for every 
n G N, the statement follows from Theorem 4 of Arcones [1]. □ 

Note that Theorem [8] can also be proven by the Central Limit Theorem for non-instanta- 
neous filters of a stationary Gaussian process given by Ho and Sun [TO] . 

We leave it as an open question whether ii < | is also a necessary condition for to be 
asymptotically normally distributed. Indeed, at least for r E Sd with d = 2 simulations 
suggest that is not asymptotically normally distributed with respect to Fh if ii > | 
(see Figure [3] below). 
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3 Estimating the probability of a change 



3.1 Ordinal patterns of order d=2 

The results of this Section are mainly based on the analysis of normal orthant probabilities 
as given by the following definition. 

Definition 9. Let n G N be fixed. For a non-singular strictly positive definite and 
symmetric matrix S G W^^^ let ■) denote the Lebesgue density of the ra-dimensional 
normal distribution with zero means and the covariance matrix S, that is 

</)(S,x) = ((27rr|S|)-^ exp{- Vs-^x} 



for X G M". We call 



<|.(S) := / 0(E,x)rfx, (li: 



[0,00 [' 

the n- dimensional normal orthant probability with respect to S. □ 

Clearly, if {Zi, Z2, . . . , Zn) is a non-degenerate mean-zero Gaussian random vector on a 
probability space {fl',A',F), then 

$((Cov(Z,,Z,))^^.=i) = P(Zi>0,Z2>0,...,Z„>0). 

The following result is well-known (see [T9j). 

Lemma 10. // {Zi, Z2, . . . , Zn) is a non-degenerate mean-zero Gaussian random vector 
on a probability space {Q',A',F) such that Va.T{Zk) > for all k eN, then 

P(Zi > 0, Z2 > 0) = 7 + 7^ arcsin pu, 

P(Zi > 0, Z2 > 0, Z3 > 0) = - + — arcsin pi2 + — arcsin pi3 

o An An 

1 



--— arcsmp23 
An 



where pij = Corr(Zj, Zj) for i,j G {1, 2, 3}. 



Note that, in general, no closed-form expressions are available for normal orthant proba- 
bilities of dimension n > A. 

Same as in Subsection 12.21 let {Xk)k€No be an equidistant discretization of fBm with the 
increment process (Ia:)^^^ = i^k — Xk-i)km- The following statement is due to Bandt 
and Shiha [3] . We refer to parts of the proof below, and thus include it here. 

Corollary 11. For He ]0, 1] one has 

iarcsin2^-i if r G {(0, 1, 2), (2, 1, 0)} 

vAH) = { 

1_ J. arcsin 2^-1 if r G {(1, 0, 2), (1, 2, 0), (0, 2, 1), (2, 0, 1)} 
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Proof. We first show the statement for r = (0,1,2). According to ([S]), if = 1 then 
P{o,i,2){H) = i. By formula ([7]) one has Varjy(lfc) = 1 for each H G ]0, 1] and /c G No, and 
hence hence CoYniXi, Yj) = CoriHiYi, Yj) for each H G ]0, 1] and i,j G N. Consequently, 
by Lemma [TO] one obtains 

= P//(ri >0,r2>0) = l + ^aTcsmpH{l) (12) 

4 /TT 

for each He ]0, 1[ with p_f/(l) as given in ([7]). Since arcsin p = 2 arcsin ^(1 + p)/2 — | 
for each p G [—1, 1], one gets p(o,i,2)(-f^) = ^ arcsin 2^~^ Now, by Lemma [3] it holds 

P(o,i,2) (H) = P(2,i,o) (H) and P(i,o,2) (^^) = P(i,2,o) (^) = P(o,2,i) (H) = P(2,o,i) (^^) 
for each H G ]0, 1], and hence the statement follows. □ 

According to Corollary [TT], the probabilities of ordinal patterns of order d = 2 are all 
monotonically related to the Hurst parameter. In particular, consider the indicator for a 
change between up to down at time k defined by 

Ck ■= l{Xfc>Xfe+i<Xfc+2} + '^{Xk<Xk+i>Xk+2} 

for G No and the probability of a change c(-) given by c{H) := Fh{Co = 1) for H G ]0, 1]. 
Because Ck = l{n2(fc)={i,o,2)} + l{n2(fc)=(i,2,o)} + l{n2(fc)=(o,2,i)} + l{n2(fc)={2,o,i)} for each 
A; G No, one has 

ciH) = P(i,o,2) (H) + P{l,2fl) (H) + P(0,2,l) (H) + P(2,0,l) (H) 

= 1 - -arcsin 2^-1 (13) 

TT 

for H G ]0, 1]. Hence, the probability of a change and the Hurst parameter H are mono- 
tonically related: the larger H, the smaller c{H). In particular, c{H) tends to | as if 
tends to 0, and c{H) = for if = 1. 

Now, consider the sample relative frequency of changes Cn given by := ^ ^^Iq for 
n G N. Since obviously 

_ ^ r W-Cn) if rG {(0,1, 2), (2, 1,0)} 

^""^ 1 \cn if rG {(1,0, 2), (1,2,0), (0,2,1), (2, 0,1)} ' 

in the case d = 2 any reasonable estimator of ordinal pattern probabilities as provided 
by Theorem [1] is an affine functional of c~„ and hence essentially has the same statistical 
properties as c~„. In particular, according to Lemma \5\ and Theorem [S], respectively, 
is a strongly consistent estimator of c(-), and asymptotically normally distributed with 
respect to for H < |. 

In the rest of this Section, we investigate the variance of c„. In Sec. [4] we consider the 
estimator of H obtained by plugging the estimate of c(-) into the monotonic functional 
relation f[T3]) . 

Notice that. Lemma [10] yields closed-form expressions also for probabilities of ordinal 
patterns of order d = 3 (see [3j). According to Theorem [H one obtains eight different 
reasonable estimators of ordinal pattern probabilites in this case where not all probabilities 
seem to be monotonically related to the Hurst parameter. 
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3.2 Variance of the sample relative frequency of changes 

Let 7h denote the autocovariance function of the stochastic process {Ck)keNo with respect 
to ¥h for H G ]0,1], that is, -fnik) = CovniCo^Ck) for H E ]0, 1] and k E Nq. Then 
Varj^(crt) can be written as 

^ n— 1 

Var^^(c„) = — (n7H(0) + 2^(n-A;)7^^(A;)). (15) 

k=l 

for 77, G N. Clearly, 7^(0) = c{H){l — c{H)), and, according to Lemma [TOl one has 

7//(i) = 2P^,(ri>o,-F2>o,r3>o)-c(i/)^ 

1 1 2 

= — arcsinpj^(2) -( arcsinpH(l)) (16) 

for H e]0,1] with pnik) as given in (JTj). For A; G N \ {1} one obtains 

^H{k) = COYH{1-Co,l-Ck) 

= FniCo = 0,Ck = 0)- FniCo = 0)^ 

= 2 FniYi >0,Y2> 0, n+i > 0, Yk+2 > 0) 

+ 2 FniY, > 0, ^2 > 0, -n+i > 0, -1^+2 > 0) 

- 4Ph(11 > 0,^2 > 0)2 (17) 

for H E ]0, 1] where, in general, no closed-form expressions are known for the probabilities 
FniY, >0,Y2> 0,n+i > 0,n+2 > 0) and Fh{Yi > 0, F2 > 0,-^+1 > 0,-^+2 > 0), 
respectively. 

Obviously, if if = 1 then c„ = P//-a.s. and hence Varj|/(c~„) = for all n G N (compare 
to ffn|) and ([S])). Note that, indeed, •ynik) = for each G Nq in this case which can be 
seen from the fact that Yi = Yj Fh-s^-s. for all i, j E N if H = 1. 

In the case H = ^ one obtains 7_h-(0) = j and •ynik) = for all k eN, and hence, according 
to formula ffTSj) . Var//(c„) = ^ for n G N. In particular, pnik) = for all A; G N if if = | 
(see ([7])), hence the Fi, F2, ^fc+i? hi (fT7|) are stochastically independent and thus 

Pj/(ri > 0, 1^2 > 0, n+i > 0, Yk+2 > 0) 

= FniY, >0,Y2> 0, -n+i > 0, -Yk+2 > 0) = PH(ri > 0, F2 > 0)^ = ^. 



Numerical evaluation of the covariances. Next, we provide a way for the numerical 
evaluation of 'ynik) in the non-trivial case H ^ {i, 1} and fc G N \ {1}. Let 7^ denote the 
set of r = (ri, r2, r3, r4) G [—1, 1]"^ such that the matrix 



S(r) 



/ 1 

ri 

r2 
\ r3 



1 

r2 



r2 

1 

ri 



^^3 

r2 
ri 
1 



is strictly positive definite. Note that, if r, s G TZ then ^Tj{h ■ r + (1 — /i) • s) x 



/ix^E(r)x+(l 
Define 



/i) x-^S(s) X > for each x G \ {0} and h E [0, 1], hence TZ is convex. 



p{{k,H,h)) := {pH{l),h-pH{k),h-pH{k + l),h-pH{k-l)) 



(19) 
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for G N \ {1}, H e ]0,1[ and h e [-1,1]. Obviously, E{p{{k, H,l))) is the co- 
variance matrix of (Yi, 1^2, ^fc+i, and S(p((fc, if, — 1))) is the covariance matrix 
of (Yi, F2, -^fc+i, -yfc+2), both with respect to Fh, and since (Fi, F2, ^fc+i, ^fc+2) and 
(Yi, Y2, — Yfc+i, — Yfc+2) are non-degenerate Gaussian with respect to if H E ]0,1[, one 
has p((fc, H, 1)) e 7^ and p{{k, H, -1)) G 7^ for each A; G N \ {1} and H G ]0, 1[. From 
the convexity of 71 it follows that p{{k, H,h)) E TZ for each k G N\{1}, H G ]0, 1[ and 
h E [—1, 1], consequently, 

v{{k,H,h)) := (<l>oS)(p((A;,if,/i))) 

is well-defined for all A; G N \ {1}, H E ]0, 1[ and h E [—1, 1] (compare to Definition [9]). 
Note that 

v{{k,H,i)) = PH(i"i>o,r2>o,n+i>o,n+2>o), 
viik,H,-i)) = p^,(ri>o,r2>o,-n+i>o,-n+2>o), 

viik,H,0)) = P^,(Fi>0,F2>0)2 

for all A; G N \ {1} and He ]0, 1[. Hence, by inserting the expressions on the left side 
into (JTTj) one obtains 



inik) 



2 ( |^((fc, ii, x)) - ^((fc, ii, -x)) ) dx (20) 



for G N \ {1} and H E ]0, 1[. The partial derivative of v with respect to h is given by 

^iik.H.x)) = pH(A:)^^|^(p((A;,ii,x))) +p^(A: + l)^^|^(p((A;,iJ,x))) 

+ p^(k-l)^^^^{piik,H,x))) 

for G N \ {1}, H E ]0,1[ and x E [—1,1]. By the reduction formula for normal orthant 
probabilities given by Plackett (see [H]), one gets the first partial derivatives of ($ o S) 
with respect to r2,r3,r4, namely 

9($oS),, 1 /II . |S(s)i3| 
s = r arcsm' 



dr2 ;,(i_, 2)1^4 2n (|S(s)n||S(s)22|)^ 

5($oS), , 1 /I 1 |S(s)23| 

arcsin 



«9r3 2n{l-sl)-2 V4 27r |S(s)22| 
^'*°^)(s) = ^fl + ia.csi„E£M| (21) 



' ' 27rn - .^2^ V4 2n |S(s 



^11 



for s = (si, S2, S3, S4) G 7?., where S(s)ij denotes the matrix obtained from S(s) by deleting 
the z-th row and j-th column of S(s). 

Formula (l20i) together with formula ( fT5l) allows to compute numerical values of Var//(c~„) 
to any desired precision. See [5] for details on the evaluation of the integral in (!20|) . The 
following Lemma will be needed below. 

Lemma 12. For each A; G No the mapping H 1-^ lH{k) is continuous on ]0, 1[. 
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Proof. For k = the statement is valid since H h-> c{H){1 — c{H)) is continuous on ]0, 1]. 
Note that H Pnik) is continuous on ]0, 1] for each A; G N (see ([7])). Consequently, 
the statement is valid for k = 1 (compare to (fT6l) ) and, furthermore, the mapping H ^ 
v{{k, H, h)) is continuous on ]0, 1[ for each G N \ {1} and h G [—1, 1]. Since ■Jnik) = 
2v{{k, H, 1)) + 2v{{k, H, -1)) - 4f ((A;, H, 0)), the statement follows. □ 

3.3 Limit behaviour 

The next Theorem establishes an asymptotically equivalent expression for 'ynik). 
Theorem 13. For each H G ]0, 1] one has 

Mk) - 4^-^if^(2if-l)2A:^^- 
7r2(l + p^(l)) 

as k tends to oo. 

Proof. First note, for H = \ and H = 1 the statement is true since both sides evaluate to 
for each A; G N (in particular Ph{^) = 1 for H = 1). Let the mapping ($ o S) : TZ ^ [0, 1] 
be given as in (fTTj) and (fT8l) . respectively. For each s = {81,82,83,34) G 71 Taylor's 
Theorem asserts the existence of some h G [0, 1] such that 

(<|.oS)(s) = ($oS)((si,0,0,0)) + X^3.^^|^((si, 0,0,0)) 

i=2 * 



2 dridr. 

i,j=2 ■< 

1 A a3($oS),, . . . ^^ 

Consequently, for each (si, S2, S3, S4) G 7?. with (si, — S2, — S3, — S4) G 7^ there exist hi G 
[0, 1] and h2 G [—1,0] such that 

($ O S)((si,S2,S3,S4)) + ($ O S)((si, -S2, -S3, -S4)) = 

2($oS)((si,0,0,0)) + ^s,s, ^ ((si, 0,0,0)) 



1 A 93($oS),, . . . ^^ 

+ o SiSjSi— — — ((si,/li-S2,/ll-S3,/ll-S4)) 



3 ^ dvidr^dri 

t,j,l=2 ■> 

1 ^ 93($oE) 

- 3 SiSjSi^^;-^^^((Si,/l2-S2,/l2-S3,/l2-S4)) . (22) 

Now, according to equation ( ITTl) . for if G ]0, 1[ and /c G N \ {1} one has 

^H{k) = 2{<^oT.){{pH{l),PH{k),pH{k + l),pH{k-l))) 

+ 2 ($ o S) ((p^,(l), -pH{k), -Pnik + 1), -ph(A: - 1))) 

-4(<l>oS)((p^(l), 0,0,0)). (23) 
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Let pnik) := H(2H - for H G ]0, 1] and A; G N \ {1}. By formula ©, pnik) is 

asymptotically equivalent to pnik) for each H G ]0, 1], hence for all i,j, I G N one has the 
asymptotic equivalence 

{pn{k)y{PH{k + i)y{PH{k-i)y ~ {pH{k)y{PH{k)y{PH{k)y 

for each H G ]0, 1]. Consequently, according to formula (l22l) and (!23|) . respectively, for 
each H e]0,l[ and /c G N \ {1} there exist hi G [0, 1] and h2 G [-1, 0] such that 

7h(A:) ~ 2(p^(fc))' ^^^^^((Ph(1),0,0,0)) + 2(pH(fc))'i?/.(^^ (24) 

where 

i,j,l=2 ■' i,j,l=2 ■' 

with p((/c, if, h)) as defined in f|T9|) . Now, by the formulas for the first partial derivatives 
of ($ o E) given in fl2T|) one obtains 

S = — r- TTT, — :^ S = — S — 



d'^r2 2ti'^{1 — s\y dr2dr^ (9r2(9r4 271^(1 — 5^ 

d^i^oT),, a2($oS),, 52($oS),, 1 

s) = — ^— ^ s — ' — 



2\ ' 



for s = (si, 0,0,0) G 7^. In particular, one has |S(s)i3| = |E(s)23| = |S(s)i4| = for 
s = (si, 0,0,0) G 7^. Putting these terms with (si, 0,0,0) = (piy(l), 0, 0, 0) into (El]), one 
gets 

7^^) - 2(p,(.))^i^^|Wj^^ 

_ 2(l-p^(l)) H\2H-lfk^''-' + 0{k^''-^RH{k)) 



vr2(l + PH(l)) 



for each H G ]0, 1[. The final claim is k^^-^ Rnik) = o{k^^-^) for each H G ]0, 1[. Let 
He ]0, 1[ be fixed. Since p((/c, H, 1)) and p((/c, H, -1)) both converge to (p//(l), 0, 0, 0) 
as k tends to oo, the convex hull of the subset 

Kh := [ [j {p((A;,ii,l)),p((A;,ii,-l))})u{(p^,(l), 0,0,0)} 

fceN\{2} 

of 7^ is closed in [-1,1]^ Note that {j^^^{p{{k,H,x)) : x G [-1,1]} is a subset of the 

a3($oS) 

dridrjdri 



convex hull of TZh- Consequently, since f .^^"f'' is continuous on IZ for all z, j, I G {2, 3, 4} 



(see [Is]), one has 

sup Rnik) < oo 

fcGN\{2} 

which shows the claim. □ 
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In the case if > | , equation (llSp together with Theorem [T2] yields simple closed-form 
asymptotically equivalent expressions for VaiHicn), namely 

Var^(c„) - ^ (25) 

where dn = -^j^+f^^-^"^ i'^-^ ~ ■ Var//(c„) tends to with a rate 

slower than in this case. Note that, from the practical viewpoint, the expressions in 
( I25l) are useful approximations of Varj/(c„) only if n is excessively large. We provide a 
way for getting more accurate practical approximations at the end of Sec. HI 

If if < I then, by equation (IT^ and Theorem [T^ one has Var//(c„) = 0(n^^). Let a 
family of mappings (/n)nGN from ]0, |[ onto M be defined by 

/„(ii) := n-VarH(c„) (26) 

for ii e ]0, 1 [ and n e N. Further, define fooiii) ■= lim^^oo fniH) for He ]0, |[. The 
following Lemma will be needed in the next section. 



Lemma 14. uniformly converges to foo on every compact subset I C ]0, |[. Moreover, 

1 4L- 



foo is continuous and strictly positive on ]0 ^ 



Proof. First, by Theorem [T3] one has Yl^=iln{k) < oo for each He ]0, |[. Hence, 
according to equation f[T^ and the dominated convergence Theorem one obtains 

n— 1 , oo ^ 

/oo(i^)-7/f(0) = lim 2 V 7^,(A;) = lim 2 V max {O, 7/f(A;) 



n— >oo * 77, n— ►OO ' ' n 

k=l k=l 

oo ^ oo 

2 V lim max{0,^^}7H(fc) = 2V7H(fc) (27) 

k=l k=l 



and, consequently, 

oo n— 1 , 



n 

k=l k=l 



oo 



2 ^ min { k/n, 1 } 7// ( A;) (28) 



fc=i 



for each H E ]0,|[. Now, let i C ]0,|[ be compact. By Lemma [T^ the mapping 
H luik) is continuous on i for every A; G N, hence one can define the family {Hk)kfz^ 
of numbers in i to be such that 'yH^ik) = Taa.x{'~fH{k) '■ H E 1} for every k E N. Then 
one has 

oo 

\foo{H) - fniH)\ < 2^min{fc/n,l}|7H,(fc)| (29) 

k=l 

uniformly for all H E I. To prove the uniform convergence of /„ to foo on i, it is sufficient 
to show that the right side of fl29l) tends to as n tends to oo. Let 5 E and define 
if* := max (i U {i + 5}) . By Theorem [13] one has \'^Hk{k)\ = 0(|7//^(/c)|), and hence 

oo CO 

^min {/c/n, 1} |7//,(A;)| = O ^ min {/c/n, l} |7/^, (A;) | j . 

k=l k=l 
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With the same arguments as in fl27p one obtains 



hm N min j/c/n, l| |7/^^ (A;) I = > hm min {A;/n, l| |7/^_^(A;)| = 

k=l k=l 

showing that the right side of (129!) tends to as n tends to oo. Now, the continuity of 
/oo follows since /„ is continuous on ]0, |[ for every n G N which is a direct consequence 
of formula f[T^ together with Lemma WI\ To establish the strict positiveness of /oo, first 
note that foo{\) = j (see Subsec. 13.21) . Further, let H G ]0, |[\{^} be fixed. According 
to Theorem [T3l one has asymptotically •ynik) > 0, hence, by equation fl28l) . there exists 
an no G N such that foo{H) - fno{H) > 0. Obviously, fn{H) > for each if G ]0, |[ and 
n G N which completes the proof. □ 



4 Confidence intervals for the ZC estimator 

According to equation (fT3!) . the probability of a change c{H) and the Hurst parameter H 
are monotonically related. Here we consider the estimator Hn of H obtained by plugging 
the estimate 'Cn of c(-) into the monotonic functional relation ( |T3|) . In order to get finite 
non-negative estimates, we define 

r log2 sin(7r(l - x)/2) + 1 for x G [O, | [ 
[ for x G [|, Ij 

and set Hn := g{cn)- Obviously, the first derivative g' of g exists and is non-zero on ]0, |[, 
namely one has 

a'U) - ^ cos(7r(l-x)/2) 

^^''^ " 21n2sin(7r(l-x)/2) ^ ^ ^''^^ 
for a; G ]0, |[. Hence, by Theorem [8] together with Theorem 2.5.2 in [16] one obtains 

{g\c{H))y\Y^iH{cn)y^HHn-H) ^ iV(0, 1) 



for each ii G ]0, |[. Note that g is continuous on [0, 1], and thus lim„^oo Hn = H P/^-a.s. 
for each H G ]0, 1] (see LemmaE]). Let fn{H) be defined as in fl26l) . According to Lemma 
[m and because the mapping H i— >• g\c{H)) is continuous on ]0, 1[, one obtains 

j.^ {9'{c{Hn))Y UHn) ^ {g'{c{H))Y U{H) ^ ^ 
™ {g'icmY fniH) {g'{c{H))f f^{H) 

P/f-a.s. for each H G ]0, |[. By Theorem 2.3.3 in [16j it follows that 

{g'{c{Hn))y\{fn{Hn)Y'^ ^{Hn- H) ^ N{Q,l) (32) 



for each H G ]0, |[. Now, let the family of mappings (s„)ngN from [0, 1] on M be defined 
by Sn{l) := 0, 

Sn{H) := n-\g'{c{H)))\^iH{cn) 
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for H G ]0, 1[, and s„(0) := ImiH^o Sn{H). Note that s„(0) is well-defined for each 
n G N since, according to Lemma [T2] and equation (pTl) . the mappings H Var//(c„) 
and H i— >■ g'{c{H)), respectively, are continuous on ]0, 1[. For n G N consider the random 
(confidence) interval 

Kn ■■= [i^„-1.96(s„(]f„))^ ]f„ + 1.96(s„(ifO)^] n [0,1]. (33) 

Let $ denote the cumulative distribution function of A^(0, 1) here. According to fl52]) one 
obtains 

¥H{HeKn) — > 2<1>(1.96)- 1 = 0.9500... (34) 

as n tends to oo. We compare this to simulations of fH{H G K^) below (see Table 

The computations of s„(if„) were carried out using formula (1371) given below. See Figure 
12] (a) for a visualization of the confidence intervals. 

Asymptotic expectation and variance. Notice that, simulations suggest that the 
actual expectation Ej:/(if„) and the actual variance Varjy (if„) of Hn asymptotically behave 
like the expectation and the variance of the asymptotic distribution of if„ obtained by 
the Delta method (see tl6j), namely, 

EniHn) ~ H-^g"{c{H))VaTH{cn) (35) 

as well as 

VarH(^n) ~ g'{c{H)fV^TH{cn) (36) 

for each H G ]0, 1] (see Table I2] below). In order to avoid confusion with the actual 
expectation and the actual variance, we call the expressions on the right side of (l35l) and 
( l36i) the asymptotic expectation and the asymptotic variance, respectively, of Hn. See 
Figure E] (b), (c) for plots of the asymptotic bias, which is the difference between the 
asymptotic expectation and H, and of the variance. The computation of Varj:/(c~„) was 
carried out using formula (137|) given below. 



(a) (b) (c) 




Figure 2: (a) Confidence intervals for the Hurst parameter obtained from Hn, for n = 128 
(black), n = 1024 (dark grey) and n = 8192 (light grey). For the value x on the abscissa, 
the ordinate shows the confidence interval Kn obtained if Hn{oj) = x. (b) Asymptotic 
bias of Hn, plotted against H. (c) Asymptotic variance of Hn, plotted against H. 
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Practical approximations. For the practical purpose of computing approximate nu- 
merical values of Var//(cn) for given G ]0, 1[ and n G N, first note that the approximation 
of 7H(fc) by the asymptotically equivalent expression 

obtained in Theorem [13] can be improved in the following way. For (ri, G such that 
(ri, r2, r2, G 7?. (with 7?. as given in the preceding section), let 

F{{ri,r2)) := $ o S((ri, r2, r2, r2)). 

Further, define 

- 4^^|i-((p„(l).0))(ff(2i/-l)i--)^' 

for i7 G ]0, 1[, m G N and > 2. Note that 7jJ^(A;) = 7H(fc) for all H G ]0, 1[ and k > 2. 
For m = 3, e = 0.01, 0.001 and several values of H, Table [T] provides the number 

/M^N • f/^o l7lr^(^) -7h(A;)I X 

that is the least k > 2 such that the relative error caused by approximating 7//(fc) by 
lip\k) is less than e. Because 7_ff(fc) and 7j^'*(fc) are asymptotically equivalent for all 
m G N, the relative error tends to as tends to oo (see [IS]), Lemma 1.1.1), that is, 
kH^\^) < oo for every e > 0. Numerical experiments suggest that J^^\k) < 'jnik) for 
all m G N, if G ]0, 1[ and k > 2. As one can see in Table [H Th {k) already yields good 
approximations of 'Jnik) for small /c if if is not too close to 1. 

H 0.05 0.15 0.25 0.35 0.45 0.55 0.65 0.75 0.85 0.95 

k^^\OM) 18 16 14 12 11 9 7 5 5 226 
fcg^ (0.001) 55 50 44 38 32 26 21 15 13 10 040 

Table 1: The numbers k\j {e) for e = 0.01, 0.001 and several values of the Hurst parameter. 



By the formulas for the first partial derivatives of ($ o S) given in (121]) one obtains 

OF 2 (II . S2(si 
"n — ~ ~ 1 7 ^ — arcsm ■ 



dr2^ ' 7r(l-si)^ V4 27r 1 + si - 2s^ 

for every s = (si,S2) in the domain of F. The partial derivatives with respect to r2 of 
higher orders are obtained by easy calculations. Those of 4th and 6th order in (p/f(l), 0) 
are given by 

d'F 4(1-p^(1))(2 + ph(1))^ 
9%^^^^^')'°)) - vr^(l + p^(l))3 ' 

d'F 16(1-ph(1))(7 + 6ph(1) + 2ph(1)T 
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For a practical approximation of Vnic^n) we propose to sum up the precise numerical 
values 7//(A;) obtained by (!20|) for /c = 2, 3, ... up to some specified n < n — 1, and then 
to take the approximate values 7//^(fc) for = + 1, . . . , n — 1, that is 

n— 1 n— 1 

Var^,(c„) ^ -(n7H(0) + 2^(n-A;)7H(A;) + 2^(n-A;)7gH^))- (37) 

A;=l A;=fi 

For the calculations behind Figure [2] we have chosen n depending on H, namely n = 
min{fcg^ (0.01), 250, n}. 

5 Simulation studies 

For the simulation of fBm we used the algorithm of Davies and Harte [8J and the random 
number generator of Mathematica 5.0.1.0. For each of the values H = 0.55,0.65,0.75, 
0.85,0.95 and the sample lengths n = 128,1024,8192, we performed 50 000 simulations 
of fBm. We only chose values H > ^, because the increments of fBm are particularly 
interesting for the modelling of long range dependence in this case (see [1]). 

5.1 Performance of the ZC estimator 

Table [2] shows simulations of Ej:/(if„) and Var/i-(if„), given by the sample mean and the 
sample variance, respectively, of the 50 000 estimates Hn obtained for simulated samples of 
length n with the true Hurst parameter H . Additionally, Table [2] provides the asymptotic 
expectation (abbreviated as. exp.) and the asymptotic variance (abbr. as. var.) of Hn 
as given by the expressions on the right side of equations (l35ll and ( |36l) . respectively. 
As one can see, the results obtained by the simulations and by equations fl35l) and fl36|) . 
respectively, agree very well. Only for H = 0.85, 0.95 the asymptotic expressions seem to 
yield slightly larger values for Varn {Hn) than the simulations. 

Table [2] also provides simulations of ^h{H G Kn), which is the probability that the 
true value of H falls within the interval Kn obtained from formula fl33l) . Theoretically, 
fniH G Kn) tends to 0.95 as n tends to cxd if if < | (see flMl) ). As one can see, for 
H = 0.55,0.65,0.75 the simulations of ^h{H G Kn) always attain values larger than 
0.95. For H = 0.85 this is still true for n = 1024 and n = 8096, while the simulations of 
FniH G Kn) are all clearly below 0.95 for H = 0.95. 

For n = 8192 and H = 0.55, 0.75, 0.95, Figure [3] shows the standardized simulated dis- 
tribution of Hn with respect to compared to the density of the standard normal 
distribution. In accordance with (|32l) . for H = 0.55 the standardized distribution of Hn 
is very close to the standard normal distribution. For H = 0.75 there is a slight difference 
between both curves, in particular, the standardized distribution of Hn seems to be more 
concentrated around than the standard normal distribution. For H = 0.95, Figure [3] 
shows a remarkable contrast between the standardized distribution of if„ and the stan- 
dard normal distribution, suggesting that if„ is not asymptotically normally distributed 
in this case. 
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0.000141 
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0.00106 








0.95 


0.951 


0.953 


0.971 


0.884 



Table 2: Simulations of KniHn), YaxniHn) and fu^H G Kn)^ obtained each for 50 000 
simulations of fBm. The rows 'as. exp.' and 'as. var.' show the asymptotic expectation 
and variance, respectively, given by the right side of equations ( l35l) and ( l36i) . respectively. 
Note that, according to (El), \i H < \ then ¥h{H G K^) tends to 0.95 as n tends to oo. 

5.2 Comparison to the HEAF estimator 

Let (yfc)fcgN be the increment process of an equidistant discretization of fBm, and define 
Yn '■= :^ ^fc=i ^A; for n G N. The HEAF estimator (HEAF meaning 'Hurst Exponent 
from Autocorrelation Function') proposed by Kettani and Gubner [T3] puts the estimate 

^ ^ E^:i(n-F„)(n+i-F„) 

of p//(l) into the monotonic functional relation Ph{^) = 2^^~^ — 1 obtained from ([7j). In 
order to receive only finite non-negative estimates of H, let 

= ^(l + log2(l + max{-l/2,p„})) 

for n G N. Table [3] shows simulations of E,h{H^^^^) and Var//(if™^^), each obtained 
for 50 000 simulations of fBm. 

As one can see, H^^-^^ has smaller variance but larger bias than the ZC estimator. For 
n = 128 the variance of the ZC estimator is about 2 times larger than the variance of 
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Figure 3: Standardized simulated distribution of if„ for n = 8192 and (a) H = 0.55 , (b) 
H = 0.75 and (c) H = 0.95, compared to the standard normal distribution (light grey). 

Hn^^^ , for n = 1024 it is about 2 to 4 times larger, and for n = 8192 it is about 2 to 8 
times larger. Notice that the highest ratios are obtained for H = 0.95. 

The loss of efficiency of if„ relative to H^^^^ is not surprising. In particular, can be 
seen as the estimator obtained by plugging the estimate sin(7r(l/2 — c„)) of puiX) into the 
monotonic functional relation pniX) = 2^^~^ — 1 (compare to (IT^ and (15U]1 . respectively), 
and clearly sin(7r(l/2 — c„)) loses efficiency relative to the estimate p„ of pij(l). However, 
because of the slow convergence of the sample mean for increments of fBm with Hurst 
parameter H > ^ (see [3]), H^^^^ has a much larger bias than the ZC estimator, in 
particular for large values of H. 

In a concluding remark, the ZC estimator of the Hurst parameter performs comparatively 
well, although it does not regard the metric structure but only the order relations between 
values of realizations. Notice that generally, for fixed r & Sd with (i G N, the computation 
of can be carried out by asymptotically 4 ({jr) (in computational steps (see [T2]). 







H 


0.55 


0.65 


0.75 


0.85 


0.95 


n = 


128 




0.538 


0.628 
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0.849 






VarH(if™^^) 


0.00377 


0.00322 


0.00267 


0.00218 


0.00167 


n = 


1024 


E^(i^™^^) 


0.548 


0.646 


0.739 


0.824 


0.893 






Var^^(^™AF) 


0.000468 


0.000399 


0.000378 


0.000373 


0.000328 


n = 


8192 


E^(i^™^^) 


0.55 


0.649 


0.746 


0.837 


0.912 






Var^,(if™AF) 


0.0000584 


0.0000515 


0.0000573 


0.0000836 


0.000101 



Table 3: Simulations of the mean and of the variance of the HEAF estimator. 
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